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Abstract 

'cK 

as 

0^ ' We investigate the mechanism of ir^i'^i (£ = e, fi) decays in which a 

^ virtual photon is emitted either from the incoming or the outgoing tt~^ . We 

point out some inconsistencies with and between two previous calculations, 
, discuss the possible experimental inputs, and estimate the branching fractions, 

f-*.. I This mechanism alone fails to explain the existing experimental data by more 

than one order-of-magnitude. But it may show itself by its interference with 
^ ■ the leading long-range mechanism dominated by the af and mesons. 

^ : PACS number(s): 12.15. y, 12.15.Ji, 13.20.Eb, 13.40.Gp 
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With the new data on the 7[^i^i~ decays coming soon [|l|, thoughts about the 

possible mechanism(s) behind them are gaining steam. One candidate is virtual photon 
emission from the incoming and outgoing Tr"*"; it was conjectured recently that its con- 
tribution to the decay amplitude may be important p. The purpose of this note is to 
scrutinize the role of this mechanism. 

The matrix element of this mechanism was given in a paper by Vainshtein, Zakharov, 
Okun, and Shifman ^ and used in the same form by Bergstrom and Singer After 
changing to contemporary notations and to the (+, — , — , — ) metrics, the matrix element 
becomes 



V2 



"'^ mj^ — ml 



2 2 



X (pk+Ptt)^ u{j)-)^^v{j)+) . (1) 

The formula we derive in this note differs from (|l|). Another formula exists in the literature 
which differs from ours "only" in the multiplicative constant. Despite the fact that our 
formula differs more substantially from Eq. (0) than it does from that in 0, we consider that 
our formula is more closely related to Eq. (|l]) than to the latter in the underlying physics, 
as we will discuss later. 

Because of these discrepancies we present the derivation of our formula in thorough 
detail. But let us first define, for completeness and later reference, the quantities which 
enter Eq. ([l|). Gp is the Fermi constant related to the weak coupling constant, and to the 
boson mass by g"^ = 4\/2Gi?m^. The decay constant of the vr"*" meson is defined by 



(0|rf(0)7'^75M(0)|p)^+ =^ F 

The definition of the decay constant, fK+, is analogous. The Vs are the elements of the 
Cabibbo-Kobayashi-Maskawa matrix. Using the and vr"*" mean lifetimes, the branching 
fractions of decays 7r+ — > /U+z/^ and K'^ — > /U+z/^ given in 0, and the method of radiative 
corrections described there, we obtained the following values: 

\f^+Vudf = (1.6418 ±0.0010) X 10-^ GeV^ 
\fK+Vus\^ = (1.2471 ± 0.0044) x 10"^ GeV^ 

which we use below in numerical calculations. The mean square radii (MSR) are defined by 
means of the low expansion of the electromagnetic form factor 

F{t) = F{0) + ^t + ..- , (2) 
o 

where F{0) = 1 (0) for charged (neutral) mesons. Finally, spinors u and v in Eq. (0) refer 
to the outgoing leptons. 

Let us now derive our formula. When the electromagnetic interaction is switched off, the 
dynamics of the system containing charged pions, kaons, and weak-gauge bosons is described 
by the effective LagrangianQ 



"^In M we used a simpler effective Lagrangian with a contact weak interaction between mesons to 
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where £free is the well-known Lagrangian of the free tt^, i^^, and fields, and 

^weak = -^^K iU^Kdd''<f>. + fK^V:,d^ct>K) + H.C. (4) 

After turning on the electromagnetic interaction by the minimal substitution — > d^+ieA'^, 
the Lagrangian acquires additional terms and becomes 

The last part contains the product of two electromagnetic field operators and does not 
operate in our case. The one-photon part contains two terms 

■^7 -^7, free ~l~ -^7, weak- 

The former was induced from the free part of the Lagrangian (^), and the latter from its 
weak interaction part @. They are given by 



P=n,K 

/:7,weak = ^A'K + + H.C. 



The amplitude of the (virtual) photon emission from W"^ is suppressed by a factor of 
{rriK/mwY , and hence, the corresponding term in C-y^free is not shown. 

To proceed further, we need to determine the interaction Hamiltonian. We can show 
that the relation 

'^int. -^weak -^7 .^77 

holds, up to non-covariant terms which vanish when the matrix element is calculated between 
the physical states P,^. After deriving the FejTiman rules in momentum space, we modified 
them by multiplying the electric charge of mesons by the corresponding form factors Fp[q^) 
[P = TT^, K~^), thereby accounting for the vr"*" and internal structure. As a result, we 
obtain the following vertex (junction) factors, starting with the well-known electromagnetic 
vertex induced from the free Lagrangian 

• —ieFp{q'^){pa + Pb)^ for the P P7 vertex. Here, Pa (Pb) is the four-momentum of 
the incoming (outgoing) meson, q = Pa — Pb is the (virtual) photon momentum, and n 
is the index connected with the photon line; 

• -ig/{2V2) Jk+V*, for the K+ W+ junction; 



show that, for the point-Uke mesons, the matrix element considered here vanishes. We checked that 
the simpler Lagrangian gives the same results as presented here, but prefer to be more rigorous. 
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• -ig/{2\/2) U+Vud for the W+ 7r+ junction; 

• ieg /{2V2) FK+{q^) fx+V*, g^"" for the K+ W+-f vertex; 

• ieg /{2V2) F^+{q^) U+Vud g^" for the W+ 7r+7 vertex. 

The contributions to the matrix element of the 7i^i~^i~ decay are depicted in Fig. 1. 

Their evaluation leads to 

U ^'^^rf T/ Vf T..^ + < FK+jt) - F^+jt) 

X (pk + p^)'' u{p-hf^v{p+) , (5) 

where t = [px — P-kY- When writing (j^) we took advantage of the fact that the contraction 
of {pk — Ptt)^ with the lepton term vanishes. Finally, if we ignore possible contributions to 
the matrix element from other mechanisms, we get the following formula for the differential 
decay rate in the mass M = \/t : 

d^K+^n+e+e- Gpa'^ |2 1 / \2 f + 

I J7T+ Vud I I Jk+ 



X \^'^{m\, ml, t) - Amj 1 + 



t 



FK+{t) - F^+{t) 



t 



(6) 



where A(x, z) = + + — 2xy — 2xz — 2yz. All quantities on the left-hand side of 
Eq. d^) are well known, except the form factors. We need their values in the time-like region 
below the physical thresholds, where they are inaccessible to direct measurement. Only for 
F.,^+{t) a small part 

Ami <t < {mx - m^Y (7) 

of the 7r+£+£~ t-interval lies in the physical region of the e^e~ annihilation. To 

determine the form factors in the region of our interest, one has to combine the experimental 
information about them in both the space-like and time-like regions with their analytic 



properties. This was undertaken, e.g., in Ref. [10|. Unfortunately, in the small-t region, 
both the TT"*" and form factors differ only a little from unity, which greatly increases 
the relative error of their difference. Therefore, one must seek alternative methods for 
determining the difference in form factors in (P); we suggest a few below. Two are based on 
the low \t\ expansion of the electromagnetic form factor (|^). Most of our methods rely on 
the vector meson dominance (VMD) description of the form factors of pions and kaons. So, 
first, we sketch the necessary VMD formulas. 

In VMD, the form factors are described by assuming that the pseudoscalar mesons couple 
to the photon via the vector meson resonances p°, a;, and (j). Let the functions V{t) {V = p, u, 
0), normalized by the condition V"(0) = 1, describe the contributions of individual resonances 
to the pseudoscalar meson form factors. The normalization conditions imposed on the form 
factors, together with the assumptions about the isospin invariance of the strong vertices, 
lead to the formulas 
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Fx+it) = Ipit) + c uj{t) + - c) m , 

FMt) = + c u{t) + - c) m , (8) 

FKo{t) = FK+{t) - F^+{t) . (9) 

In principle, the constant c can be fixed by assuming the (flavor) SU(3) invariance. Because 
the latter is violated, we will, instead, consider c a phenomenological parameter, and flx its 
value by experimental data. But before proceeding further, we must specify the functions 
V{t). For narrow resonances uj and 0, which, in addition, do not have any open decay 
channels in our t-range, we can use the prescription 

vit) = , (10) 



m 



V 



which originates from the free-vector-particle propagator, taking into account that the term 
in its numerator containing four-momenta does not contribute. 

But the situation with the p° is more complicated. This is evident from the fact that 
using the simple formula (|10D for p{t) gives the MSR of 71+ equal to (0.3940 ± 0.0008) fm^, 
which disagrees with the experimental flndings that give (0.439 ± 0.008) fm^ [|r^. We will, 
therefore, use the form that properly accounts for the dynamics of resonances (see, e.g., |1T2[ ) 



= 'f «) , (U) 

where m'^p{t) is the running mass squared, rrip is the resonant mass, and Tp{t) is the variable 
total width. The latter vanishes in our kinematic range, with the exclusion of the small 
region already mentioned (0). We take the function m'j^{t) from our recent work |jl2|; this 
gives the vr"*" MSR of (0.446 ± 0.006) fm^, in good agreement with data. 

Below, we briefly describe six methods for evaluating the form-factor difference, which 
we need to be able to calculate the branching fraction (^ of the K+ TT+i+i~ decays. The 
corresponding results for both the e+e~ and p+ pT modes are shown in Table 1. 

1. Both F^^+it) and F^+it) are taken from Ref. []TU[. We do not have any way to assess 
the errors, so the results are shown with the "approximate" sign. 

2. Formula (^ is used, together with the parametrization of FKc{t) from [0. The same 
caveat as above applies. 

3. Low-t expansion (§) is used with 

(r2)^+ - (r2)x+ = (0.100 ± 0.045) fm^ , (12) 



taken from Ref. [13 . 
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4. The formula @ and parametrization (|]) for Fxo{t) are used with 

(r2)^o = -(0.054 ± 0.026) fml (13) 

This value was determined in |T^ by measuring the coherent regeneration of KgS by 



atomic electrons; the method was proposed by Zel'dovich [|T5 



5. Formulas (^ and (^ are used with c = 0.07±0.29, which was determined by matching 
the MSR difference (1121). 



6. As above, but the MSR of (0) was matched with c = 0.36 ± 0.17. 

Of all of those methods, we believe that the last two are the most reliable. As we argued, the 
calculation based directly on the form factors (method 1) is hampered by the large relative 
errors in their difference. Method 2, using the form factor, is probably more reliable, but 
similarly here, we cannot estimate the final error. Methods 3 and 4 lead to a constant form 
factor of the — > 7r~^ decay, defined below in Eq. (0), which is a bad approximation 
T^ . Our preference for methods 5 and 6 stems from the fact that the systematic errors of 



the VMD method are certainly smaller than the errors of the data we use as input. Thus, 
the extrapolation from t = 0, where we fixed parameter c [see Eq. (^], to higher values of 
t is well defined. Morever, c determines the relative weight of the a; and contributions, 
which do not differ greatly in our region of t. So, the value of c is not critical. 

The branching fractions are not the only quantities that can be compared to experimental 
data. The decay form factor f{t), which also is very important, is defined through the 
matrix-element parametrization 

M = C fit) (pK + p.r uiP-hAP+) (14) 



and the normalization condition /(O) = 1. Comparing Eqs. (^ and (|I4D reveals that the 
mechanism considered here exhibits a form factor 

l i^x+(t)-F.+ (t) 

In Fig. 2, we plot the absolute value of this form factor assuming the VMD relations (H) and 
(^) with c = 0.36 (method 6), and of the form factor calculated for the long-range mechanism 
dominated by the af and mesons |l7|,[T2|. What is finally observed experimentally is the 
form factor belonging to the superposition of those two mechanisms. This most important 
result is depicted by the solid line. Because the normalisation of the meson-dominance 
matrix element is yet a little uncertain, see discussion in , we considered its multiplicative 
constant a free parameter. Its value was chosen to get the correct branching fraction for the 
e+e" mode (2.74 x 10"'', 0). For comparison, we also show the linear parametrization of 
the form factor used by experimentalists 

fit) = 1 + , 

where A = 0.105 [0. Taking into account the experimental errors of A, which are 



0.035 (stat.) and 0.015 (syst.) [TO], we can show that the solid curve is compatible with 
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the published experimental results. However, the preliminary data of the BNL-E865 col- 
laboration [m indicates a steeper slope, which may be a problem for the meson dominance 
model, even when it is supplemented with the mechanism considered in this note. 

Having fixed the form factor that belongs to the superposition of the two mechanisms 
mentioned above, we can evaluate the fi~^ fi~ / e'^ e~ branching ratio. It comes out as 



B{K+ TT+e+e-) 



0.248 ± 0.002. 



The quoted error reflects the uncertainties of the parameter c and of the running mass 
squared |12| in our t-region. Using the recommended value (2.74 ± 0.23) x 10~^ of the 
branching fraction of the e~^e~ mode [Q, we end up with the prediction 



n fi fi 



(6.8 ±0.6) X 10^ 



Let us compare now our formula for the matrix element and the numerical results with 
those of previous studies 

Our matrix element (^) differs from that derived in 0, see Eq. (|I[), even after the form- 
factor difference in our formula is expressed in terms of the mean radii squared, using Eq. (^. 
We could not locate the source of the discrepancy. To see how serious this discrepancy is 
from a pragmatic point of view, we calculated the absolute value of the matrix element ratio 
using the mean square radii of vr"*" and of 0.44 fm^ and 0.34 fm^, respectively [IT, 
The result is 



M (this work) 



M{Ref. i; 



mj^ + mi 



0.18 . 



A formula for the matrix element of the same "inner-bremsstrahlung" mechanism was 
proposed also in It contains a simple difference of the vr"*" and electromagnetic 
form factors, in agreement with our findings and in disagreement with ([l|). However, the 
same form was obtained by different means. The authors of simply assumed that the 
— »• 7r~^ transition amplitude was a constant, whereas in [Q, and also in our approach, 
this amplitude is proportional to the momentum squared. The latter is equal to the mass 
squared of that of the two mesons which is on the mass shell. In our approach the form factor 
difference results from the interplay between the momentum-dependent transition amplitude 
and the contributions to the matrix element that originate in the weak-interaction part of 
the Lagrangian (^). 

The multiplicative constant in the matrix element presented in is fixed, through a 
chain of reasoning, by the data on the K ^ titi and A^^^'^t decay rates. This shows 

that the physics assumed to lie behind the 7r+ transition is different from that here 

and in Ref. Q], where the former source of information is replaced by the 7r+ — > yU'^z/^ decay 
rate. 

What concerns the magnitude of the matrix element here and in 0, their estimate 



[7r+\H^\K+) 



3.9 X 10"^ GeV^ 



should be compared to ours 
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Gf 
2V2 



l/.+ Kd/x+K.! {ml + ml) = 4.91 x IQ-^ 



The strikingly different magnitudes of the matrix element lead to different judgements about 
the role of the considered mechanism. In it was deemed very important, with the matrix 
element twice as big as that required by experimental branching fraction for the e~^e~ mode. 
Here, we found (see Table 1) that this mechanism, if taken alone, leads to branching fractions 
that are at least one order-of-magnitude below the experimental values. 

It was shown in ^ that the ai/p-meson dominance makes the leading contribution to 
the matrix element of the —>■ 7T^i~^i~ decays. Now we have shown that the matrix 
element of the mechanism considered here is several times smaller (4-8 times, if we quote 
the results of methods 5 and 6, which we trust the most). Neverthless, its interference with 
the dominant mechanism somewhat improves the behavior of the form factor, which is a 



little too flat for the dominant mechanism alone [112 
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TABLES 



TABLE L Branching fractions of the — > 7r~^e~^e and n^n^fx decays calculated 

using different inputs and compared to data (1) X+ and 7r+ form factors |jl^; (2) form 
factor 0; (3) MRS difference of K+ and tt+ 0; (4) MRS of K° [|14); (5) VMD form factor of 
K° with c fixed by the difference in K+ and 7r+ MRS |3|; (6) as in (5), but fixed by the MRS 

& 



Method 




1 


2 


3 


4 


5 


6 


Exp. data 


B{e+e-) 


X 10^ 
X 10^ 


L3 
« 0.3 


PS 22 
f« 7.4 


1 q+2.1 


n ^4+0-65 
U-C»^-0.40 


14+- 

^•"-2.4 


4 6+^-^ 

^•"-2.9 


274 ± 23 
50± 10 
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FIGURES 




FIG. 1. Fcynman diagrams contributing to the matrix clement of the n^i'^i^ decay. 

Contributions (la) and (lb) are proportional to the form factor, (2a) and (2b) to the tt^ form 
factor. Contributions (la) and (2a) were generated from the free part of the Lagrangian, and (lb) 
and (2b) from its weak-interaction part. Contribution (3) is suppressed by an additional factor of 
{rriK/mw)'^ and therefore is ignored. 
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FIG. 2. Form factor of the decay K'^ ir^i^i^: The mechanism considered in this work 
taken alone (dashed hne); The form factor coming from the ai/p dominance l?],^ (dash-dotted); 
Superposition of the oi /p dominance and the present mechanism (sohd); The hnear parametrization 
used to fit data |jl^ (dotted). 
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